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Abstract 

An algorithm for the calculation of the resistance between two arbitrary nodes in an 
arbitrary distance-regular resistor network is provided, where the calculation is based on 
stratification introduced in [1] and Stieltjes transform of the spectral distribution (Stielt- 
jes function) associated with the network. It is shown that the resistances between a 
node a and all nodes (3 belonging to the same stratum with respect to the a {R a p(i) , (3 
belonging to the i-ih stratum with respect to the a) are the same. Also, the analytical 
formulas for two-point resistances R aj3 (i) , i = 1,2,3 are given in terms of the the size of the 
network and corresponding intersection numbers. In particular, the two-point resistances 
in a strongly regular network are given in terms of the its parameters (v, k, A,/x). More- 
over, the lower and upper bounds for two-point resistances in strongly regular networks 
are discussed. 

Keywords:two-point resistance, association scheme, distance-regular net- 
works, Stieltjes function 

PACs Index: 01.55.+b, 02.10.Yn 
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1 Introduction 

A classic problem in electric circuit theory studied by numerous authors over many years is 
the computation of the resistance between two nodes in a resistor network (see, e.g., [2]). 
Besides being a central problem in electric circuit theory, the computation of resistances is 
also relevant to a wide range of problems ranging from random walks (see [3]), the theory of 
harmonic functions [4], to lattice Greens functions [5, 6, 7, 8, 9]. The connection with these 
problems originates from the fact that electrical potentials on a grid are governed by the same 
difference equations as those occurring in the other problems. For this reason, the resistance 
problem is often studied from the point of view of solving the difference equations, which is 
most conveniently carried out for infinite networks. In the case of Greens function approach, 
for example, past efforts [2], [23] have been focused mainly on infinite lattices. Little attention 
has been paid to finite networks, even though the latter are those occurring in real life. In this 
paper, we take up this problem and present a general formulation for computing two-point 
resistances in finite networks. Particularly, we show that known results for infinite networks 
are recovered by taking the infinite-size limit. 

The study of electric networks was formulated by Kirchhoff [11] more than 150 years ago 
as an instance of a linear analysis. Our starting point is along the same line by considering the 
Laplacian matrix associated with a network. The Laplacian is a matrix whose off-diagonal en- 
tries are the conductances connecting pairs of nodes. Just as in graph theory where everything 
about a graph is described by its adjacency matrix (whose elements is 1 if two vertices are 
connected and otherwise), everything about an electric network is described by its Laplacian. 
The author of [12], has been derived an expression for the two-point resistance between two 
arbitrary nodes a and (5 of a regular network in terms of the matrix entries L~*, Lzl and 
L~p, where L~ l is the pseudo inverse of the Laplacian matrix. Here in this work, based on 
stratification introduced in [1] and spectral analysis method, we introduce an procedure for 
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calculating two-point resistances in distance-regular resistor networks in terms of the Stieltjes 
function G^{x) associated with the adjacency matrix of the network and its derivatives. Al- 
though, we discuss the case of distance-regular networks, but also the method can be used for 
any arbitrary regular network. It should be noticed that, in this way, the two-point resistances 
are calculated straightforwardly without any need to know the spectrum of the network. Also, 
it is shown that the resistances between a node a and all nodes (3 belonging to the same stra- 
tum with respect to the a (R a p(i) , (3 belonging to the i-th stratum with respect to the a) are 
the same. We give the analytical formulas for two-point resistances R a/3 (i),i = 1,2,3 in terms 
of the network's characteristics such as the size of the network and its intersection array. In 
particular, the two-point resistances in a strongly regular network are given in terms of the 
network's parameters (v, k, A, //). Moreover, we discuss the lower and upper bounds for two- 
point resistances in strongly regular networks. From the fact that, the two-point resistances on 
a network depend on the corresponding Stieltjes function G^{x) and that G^(x) is written as a 
continued fraction, the two-point resistances on an infinite-size network can be approximated 
with those of the corresponding finite-size networks. 

The organization of the paper is as follows. In section 2, we give some preliminaries such 
as association schemes, distance-regular networks, stratification of these networks and Stieltjes 
function associated with the network. In section 3, two-point resistances in distance-regular 
networks are given in terms of the Stieltjes function and its derivatives. Also, the resistances 
R a p(i), for i — 1, 2, 3 are given in terms of the network's intersection array. In particular, two- 
point resistances in a strongly regular network are given in terms of the network's parameters, 
also lower and upper bounds for the two-point resistances in these networks are discussed. Sec- 
tion 4 is devoted to calculating two-point resistances R aj3 (i) for i — 1, 2, 3 in some important 
examples of distance-regular networks, such as complete network, strongly regular networks 
(distance-regular networks with diameter 2), e.g. Petersen and normal subgroup scheme net- 
works [1], (i-cube (d dimensional hypercube) and Johnson networks. The paper is ended with a 
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brief conclusion and an appendix containing a table for two-point resistances R a p(i), % = 1, 2, 3 
of some important distance-regular resistor networks with size less than 70. 

2 Preliminaries 

In this section we give some preliminaries such as definitions related to association schemes, 
corresponding stratification, distance-regular networks and Stieltjes function associated with 
the network. 

2.1 Association schemes 

First we recall the definition of association schemes. The reader is referred to Ref.[13], for 
further information on association schemes. 

Definition 2.1 (Symmetric association schemes). Let V be a set of vertices, and let Ri(i = 
0, 1, d) be nonempty relations on V (i.e., subset of V x V). Let the following conditions (1), 
(2), (3) and (4) be satisfied. Then, the relations {Ri}o<i<d on V x V satisfying the following 
conditions 

(1) {Ri}o<i<d is a partition of V x V 

(2) Rq = {(a, a) : a G V} 

(3) Ri = R\ for < % < d, where R\ = {{(3, a) : {a, (3) G R { } 

(4) For (a,/3) G Rk, the number ■ =| {7 G X : (a, f3) G Ri and (7, /3) G | does not 
depend on (a, (5) but only on i,j and k, 

define a symmetric association scheme of class d onV which is denoted by Y = (V, {Ri}o<i<d}- 
Furthermore, if we have = p^ for all k = 0, 2, d, then Y is called commutative. 

The intersection number can be interpreted as the number of vertices which have relation 
i and j with vertices a and (3, respectively provided that (a, (3) G Rk, and it is the same for 
all element of relation R k . For all integers i (0 < i < d), set ki — p\ and note that fc, 7^ 0, 
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since Ri is non-empty. We refer to hi as the i-th valency of Y . 

Let Y = (X, {Ri}o<i<d) be a commutative symmetric association scheme of class d, then 
the matrices A , A ly A d defined by 

f 1 if (a, P) G Ri 
{A*)*j> = , (2-1) 

1 otherwise 

are adjacency matrices of Y such that 

^ = Y.V%A k . (2-2) 

From (2-2), it is seen that the adjacency matrices Aq, Ai, ...,Ad form a basis for a commutative 
algebra A known as the Bose-Mesner algebra of Y. 

Definition 2.2 (P-polynomial property) Y = (X, {Ri} 0<i<d ) is said to be P-polynomial (with 
respect to the ordering R Q , R d of the associate classes) whenever for all integers h, i,j (0 < 

h,i,j < d), 

Pij = if one of h, i,j is greater than the sum of the other two, 

p\j 7^ if one of h,i,j equals the sum of the other two. (2-3) 

It is shown in [14] that in the case of P-polynomial schemes, A; L = Pi(A) (0 < i < d), where 
Pi is a polynomial with real coefficients and degree exactly i. In particular, A multiplicatively 
generates the Bose-Mesner algebra. 

Finally the underlying graph of an association scheme T = {V, Ri) is an undirected con- 
nected graph, where the set V and Ri consist of its vertices and edges, respectively. Obviously 
replacing Ri with one of other relation such as Ri, for i ^ 0, 1 will also give us an underlying 
graph T = (V, R^ (not necessarily a connected graph ) with the same set of vertices but a new 
set of edges Ri. 
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2.2 Stratifications 

For a given vertex a G V, let Ri(a) := {j3 G V : G denote the set of all vertices 

having the relation Ri with a. Then, the vertex set V can be written as disjoint union of Ri(a) 
for % — 0, 1, 2, d, i.e., 

V=(j Ma). (2-4) 

i=0 

We fix a point o G V as an origin of a distance regular graph (the underlying graph of an 
association scheme), called reference vertex. Then, the relation (2-4) stratifies the underlying 
graph into a disjoint union of associate classes Ri(o). With each associate class Ri(o) we 
associate a unit vector in l 2 (V) defined by 

\*i) = A= £ l«>> (2-5) 

V Kt aeRi(o) 

where, \a) denotes the eigenket of a-th vertex at the associate class Ri(o) and Ki = \Ri(o)\ is 
called the i-ih valency of the graph (i.e., Ki = The closed subspace of 1 2 {V) spanned by 
is denoted by A(G). Since becomes a complete orthonormal basis of A(G), we 

often write 

A(G)=^©C|0 l ). (2-6) 

i 

Let Ai be the adjacency matrix of a graph T = (V, R) for reference state \<f> ) (\<f> ) = |o), with 
o G V as reference vertex), we have 

AM= E 1/5)- (2-7) 
Then by using unit vectors (2-5) and (2-7), we have 

Ai\<f> ) = ^l\<pi)- (2-8) 



2.3 Distance-regular graphs 

Here in this section we consider some set of important graphs called distance regular graphs, 
where the relations are based on distance function defined as follows: A finite sequence 
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ao, oti, a n G V is called a walk of length n (or of n steps) if a k -i ~ oik for all k = 1, 2, n. 
For a 7^ /3 let d(a, (3) be the length of the shortest walk connecting a and f3, therefore d(a, (3) 
gives the distance between vertices a and (3 hence it is called the distance function and we 
have d(a,a) = for all a G V and d(a,j3) = 1 if and only if a ~ (3. Therefore, the distance 
regular graphs become metric spaces with the distance function d. 

An undirected connected graph V = (V, Ri) is called distance regular graph if it is the 
underlying graph of a P-polynomial association scheme with relations defined as: G Ri 

if and only if d(a,(3) = i, for % — 0, 1, ...,d, where d :=m&x{d(a, (3) : a, (3 G V} is called the 
diameter of the graph. Usually in distance regular graphs, the relations Ri are denoted by IV 

Now, in any connected graph, for every (3 G Ri(a) we have 

Ri{(3) C i2i_i(a) U i2i(a) U R i+1 (a). (2-9) 

Hence in a distance regular graph, = (for i ^ 0, j is not — 1, i, i + 1}). The intersection 
numbers of the graph are defined as 

Oi=p\i, bi=p\ +1 ,i, c i =p\_ lil . (2-10) 
The intersection numbers (2-10) and the valencies Ki satisfy the following obvious conditions 

O'i bi + Cj = Ki_\bi_i = /tjCj, i = 1, d, 

k = ci = 1, 6 = «i = (c = b d = 0). (2-11) 

Thus all parameters of the graph can be obtained from the intersection array {b , b d _i, ci, q}. 

It should be noticed that, for distance-regular graphs, the unit vectors \4>i) for i — 0,1, d 
defined as in (2-5), satisfy the following three-term recursion relations 

A\<f>i) = A+il^+i) + aM + (2-12) 
where, the coefficients ai and (3i are defined as 

a k = K-b k -c k , uj k = Pi = &fc_ic fc , fc = l,...,tZ, (2-13) 
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i.e., in the basis of unit vectors {\4>i),i = 0, 1, ...,d}, the adjacency matrix A is projected to 
the following symmetric tridiagonal form: 

/.on n \ 



a A 
/3i cti /3 2 

#2 «3 & 










(2-14) 



... (3 d -i a d ~i (3 d 
V ... (3 d a d J 

In Ref. [20], it has be shown that, the coefficients on and $ can be also obtained easily by 
using the Lanczos iteration algorithm. Hereafter, we will refer to the parameters ol; l and Ui in 
(2-13) as QD (Quantum decomposition) parameters. 

One should notice that, in distance regular graphs stratification is reference state indepen- 
dent, namely we can choose every vertex as a reference state, while the stratification of more 
general graphs may be reference dependent. 



2.4 Stieltjes function associated with the network 

It is well known that, for any pair (A, |0 O )) of a matrix A and a vector |0 O ), it can be assigned 
a measure /i as follows 

u.(x) = (<p \E(x)\<( )0 ), (2-15) 

where E(x) = J2i \ui){ui\ is the operator of projection onto the eigenspace of A corresponding 
to eigenvalue x, i.e., 

A = J xE(x)dx. (2-16) 
It is easy to see that, for any polynomial P{A) we have 



P(A) = J P(x)E(x)dx, 



(2-17) 
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where for discrete spectrum the above integrals are replaced by summation. Therefore, using 
the relations (2-15) and (2-17), the expectation value of powers of adjacency matrix A over 
starting site \<f> ) can be written as 

(0o|^ m |0o) = / x m fi{dx), m = 0,1,2,.... (2-18) 

JR 

The existence of a spectral distribution satisfying (2-18) is a consequence of Hamburgers the- 
orem, see e.g., Shohat and Tamarkin [[15], Theorem 1.2]. 

Obviously relation (2-18) implies an isomorphism from the Hilbert space of the stratification 
onto the closed linear span of the orthogonal polynomials with respect to the measure fi. It is 
well known that [20] for distance-regular graphs we have 

|0 i ) = P i (A)|0 o >, (2-19) 

where P; is a polynomial with real coefficients and degree i. Now, substituting (2-19) in (2-12), 
we get three term recursion relations between polynomials Pj(A), which leads to the following 
three term recursion between polynomials Pj(x) 

xP k {x) = !3 k+1 P k+1 {x) + a k P k {x) + p k P k ^(x) (2-20) 

for k = 0, d - 1, with P (x) = 1. 
Multiplying by [3i.../3 k we obtain 

Pi...p k xP k {x) = fa...p k+1 P k+1 {x) + a k (3 l ...(3 k P k {x) + (3 2 k .(3 1 ...(3 k - 1 P k - 1 {x). (2-21) 

By rescaling P k as Q k = (3i...(3 k P k , the spectral distribution \i under question is characterized 
by the property of orthonormal polynomials {Q k } defined recurrently by 

Q (x) = l, Qi(x)=x, 

xQ k {x) = Q k+l {x) + a k Q k {x) + PlQk-^x), k > 1. (2-22) 
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If such a spectral distribution is unique, the spectral distribution /i is determined by the 
identity 

G f u(dy)_ = 1 = ggi(x) = ^ At 

" JRX-y z _ aQ PI Q d (x) f^x-xi' 

z ~ ai 

z — an 

z Z — OLQ 

where, x t are the roots of polynomial Qd(x). G^z) is called the Stieltjes/Hilbert transform of 
spectral distribution /j, or Stieltjes function and polynomials {Q^} are defined recurrently as 

QS\x) = l, Q*t\x)=x-a l , 

xQ£\x) = Q^Ux) + a k+l Q k 1] {x) + (3 2 k+1 Q^Ux) , k > 1, (2-24) 
respectively. The coefficients A x appearing in (2-23) are calculated as 

A x = lim(2-xi)G>0z). (2-25) 

Z *Xi 

(for more details see Refs.[15, 16, 17, 18].) 

3 Two-point resistances in regular resistor networks 

A classic problem in electric circuit theory studied by numerous authors over many years, is 
the computation of the resistance between two nodes in a resistor network (see, e.g., [2]). The 
results obtained in this section show that, there is a close connection between the techniques 
introduced in section 2 such as Hilbert space of the stratification and the Stieltjes function 
and electrical concept of resistance between two arbitrary nodes of regular networks and these 
techniques can be employed for calculating two-point resistances. 

For a given regular graph Y with n vertices and adjacency matrix A, let = rji be the 
resistance of the resistor connecting vertices % and j. Hence, the conductance is = r' 1 = 
so that Cij = if there is no resistor connecting i and j. Denote the electric potential at the 
i-th vertex by Vi and the net current flowing into the network at the i-th vertex by ij (which 
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is zero if the i-th vertex is not connected to the external world). Since there exist no sinks 
or sources of current including the external world, we have the constraint Yh=i h — 0. The 
Kirchhoff law states 

E c lJ {V l -V J )=I l , i = l,2,...,ra. (3-26) 

Explicitly, Eq.(3-26) reads 

LV = I, (3-27) 
where, V and / are n- vectors whose components are Vi and I iy respectively and 

^ = E c il0(<l-E c «l < >01 (3-28) 



is the Laplacian of the graph T with 



Ci= E c ^ ( 3 - 29 ) 



for each vertex a. Hereafter, we will assume that all nonzero resistances are equal to r, then 
the off-diagonal elements of — L are precisely those of ^A, i.e., 

L — -(kI - A), (3-30) 

r 

with k = deg(a), for each vertex a. It should be noticed that, L has eigenvector (1, 1, 1)* 
with eigenvalue 0. Therefore, L is not invertible and so we define the psudo-inverse of L as 

L- 1 = £ A- 1 ^, (3-31) 

where, E { is the operator of projection onto the eigenspace of L^ 1 corresponding to eigenvalue 
Aj. Following the result of [12] and that L^ 1 is a real matrix, the resistance between vertices 
a and j3 is given by 

R af3 = (alL'^a) - 2(a\L- 1 \{3) + {^L' 1 ^). (3-32) 

In this paper, we will consider distance-regular graphs as resistor networks. Then, the diag- 
onal entries of Lr x are independent of the vertex, i.e., = for all a, (3 G V. Therefore, 
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from the relation (3-32), we can obtain the two-point resistance between two arbitrary nodes 
a and [3 as follows is written as 

R a p = 2{L-l-L-}). (3-33) 

Now, let a and f3 be two arbitrary nodes of the network such that (3 belongs to the /-th 
stratum with respect to a, i.e., (5 G Ri(a) (we choose one of the nodes, here a, as reference 
node). Then, for calculating the matrix entries and Lp* in (3-32), we use the Stieltjes 
function to obtain 

£ - = ^Tih^ = 'U, =T %^ = r ^ (GM ~ (3 - 34) 



and 



= rip — -a) = — — (07 — -a) = — — (a — 



a) 



T f W ^)~E W . (3-35) 



■/#-{/*} K-X y/Ei K - X. 

where, we have considered xo = k (k is the eigenvalue corresponding to the idempotent E ). 
Then, by using (3-33), the two-point resistance R a p(i) in the network is given by 

For evaluating the term Y,i,i^o Ai ^-x^ in (3-36), we need to calculate 



A x m 

I m :=J2^Z^i for rn = 0,1,..., I. (3-37) 

i,i^0 K x * 



To do so, we write the term (3-37) as 



l = Axi = Aidxi - n) m - E^-iycr^r') 



that is, we have 



- £ - K ) m - 1 - E(-i)^N E (3-38) 

i,i^0 1=1 i,i^0 K x% 



m—l m 

/m = -E (-i)'cr E 4*r ,_1 - E(- l )' c 'N^- (3-39) 

Z=0 i,i^0 i=l 
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Therefore, I m can be calculated recursively, if we able to calculate the term J2i,#o M^T" 1-1 
for I = 0, 1, ...,m — 1 appearing in (3-38). For example, for m — 1, we obtain 

h = E = -^Ai + 4 = -l+4 + ^ (3-40) 

In order to evaluation of the sum Y,i,i^o Aixf , we rescale the roots $i clS 1 where £ is some 
nonzero constant. Then, we will have 

Now, we take the m-th derivative of (3-41) to obtain 

where, at the limit of the large x, one can obtain the following simple form 



2m ^s\g,mo) = m!( g^g±^ ) ,3-43) 



Therefore, we obtain 



^ A ^ = ^ Jl^ m+ " ( 3 " 44 ) 
3.1 Two-point resistances up to the third stratum 

In this subsection we give analytical formulas for two-point resistances R a/3 (i),i = 1,2,3, in 
terms of intersection numbers. 

It should be noticed that, for two arbitrary nodes a and (5 such that (5 G Ri(a), we have 
P±(x) = Therefore, by using (3-35) and (3-40) we obtain 

E^ = 7E^. + 'E^ (3-45) 

Therefore, by using (3-33), we obtain the following simple result for all (3 G Ri(a) 

2r — 2r 2r 1 

RafW = ~ E A = -(i - A)) = -(i - -), (3-46) 

re i ^ re rev 
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where, v is the number of vertices of the graph, and in the last equality we have used the fact 
that for regular graphs, we have 

A = -. (3-47) 

v 

In the following, we give analytical formulas for calculating two-point resistances R a s(2) and 
R a p(3), where R a p(2) (R a p( 3 ) ) denotes the mutual resistances between a and all j3 G R2{ot) (all 
(5 G R 3 (a)). 

By using (2-22) and that P k = ^== , we have P 2 (x) = -j==(x 2 - a ± x -u)\). Then, from 
(3-35) after some simplifications we obtain for (3 G R2(o) 



v A ■ 

L tpm = /rn-nr ^- E A ^ + ("i - K ) E A i + < K ~ «i - x ) E — ^7 

By substituting «i = k — bi — c\ in k(k — a± — 1), we obtain 



k(k — «! — !) = k(6i + ci — 1) = Kb\. (3-49) 



Then, the coefficient of the term J2i,i^o i n (3-48) is 



rKh rKh =J^ = r. (3-50) 



Therefore, (3-48) can be written as 

L-; (2) = -j=f={~ E + («i - «) E A) + r £ (3-51) 
where, the sum X)i,i^o ^i 2 ^ can be calculated by using (3-44). It can be easily shown that 

J™ lx 2 ^(\G,(x/0)} = a^ 2 - 6d-i, (3-52) 



where, a d _ 2 and are the coefficients of x d ~ 2 and a^" 1 in an d respectively. From 
the recursion relations (2-22) and (2-24), one can see that ad-2 = bd-i = —{ot\ + ... + ad). 
Therefore, from (3-44) and (3-52) we obtain 



A i x i = ~A K = --. (3-53) 
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Then, by using (3-32) and (3-51), one can write R a p(2) as follows 

_. 2r 2k — eii 

R a/3 (2) = {(« - an) }, (3-54) 

y /UJ 1 UJ2K2 V 

where, by using (2-11) and (2-13), we obtain the following main result in terms of the inter- 
section numbers of the graph 

^ Q /3(2) = + 1 1- (3-55) 

Now, consider /3 e -^3(0;). Then, by using (2-22) and P k = p® k p k we obtain P 3 (x) = 
-j=== [x 3 — («i + ol2)x 1 — [uj\ + oj<i — 01x012)% + «2^i)- As above, after some calculations, we 
obtain for (3 G R 3 (a) 

2 

L~l 3) = { 2(a d _ 3 - b d _ 2 + - 6 d _ia ci _2) - («i + a 2 - k) — 

a ^ ^^2^3 K 3 V V 

V — 1 A 

(/t 2 -/t(«i+o;2)-^i-^2+aia2)( ) + (/t 3 -/t 2 (o;i+a2)-^(^i+^2-ttia2)+tt2^i) X! " — }■ 

V i,i^o k- Xi 

(3-56) 

Again, by substituting eti, a 2 , o;i and c^2 from (2-13), we have 

(k 3 — K 2 (ai + a 2 ) — + ^2 — ol\ol 2 ) + «2^i) = , = = 1. (3-57) 



y/UJiUJ20JzKz ' V r Kb 1 ~C2~b 2 ~C^K 3 

Therefore, (3-56) can be written as follows 

2 

L~l 3) = { 2(a d _ 3 - b d _ 2 + b\_ x - b^a^) - («i + a 2 - «) 

a P y/u}iU 2 UJ 3 K 3 V V 

v — 1 .---> A- 

(k 2 - k{ol\ + a 2 ) -UJ1-UJ2 + ot 1 a 2 )( )} + r ^ ~- (3-58) 



In (3-56), we have used the following equality 



lim [ x s ^(Ig^x/O)] = 2(a d _ 3 - &d-2 + - ^-1^-2), (3-59) 

where, a^_ 3 and 6 d _ 2 are the coefficients of x d ~ 3 and x d ~ 2 in an d respectively. From 
the recursion relations (2-22) and (2-24), one can see that a d _ 3 = llf<j=i otiOtj — (uo\ + ...+Ud-i) 
and b d _ 2 = a d ^ 3 - u d . Therefore, we have a d _ 3 - b d _ 2 + _ & d-i«d-2 = ^d- 
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Then, by using (3-32), R aj 3(3) is given by 

R a0 (3) = -{uj d + (a 1 + a 2 -2K,)- + (K 2 - K,(a 1 + a 2 ) -mi -^2 + ^10:2) ( )}• 

In terms of the intersection numbers of the graph, we obtain the following main result 



7? ^ Ih r h+r-Mhh (^0 + l)fe + C 2 )+&l(&0 + 6 2 ) 1 

-^a/3 (a) = mri t>d ~ lCd + 02-00 + C2 + bib 2 }. (3-6i) 

3.2 Two-point resistances in infinite regular networks 

As the results (3-33) and (3-44) show, the two-point resistances on a network depend only 
on the Stieltjes function G^{x) corresponding to the network. Clearly, each Stieltjes function 
corresponding to an infinite network has a unique representation as an infinite continued 
fraction as follows 

GM = / ^ = ^ p , (3-62) 

Jrx-v z _ ao ^ 

— 



z-a 3 

where, the sequence ceo, &i, Pi, P2, ■■■ never stops. It is well known that all infinite continued 
fraction expansions as in (3-62) can be approximated with finite ones. Therefore, the two- 
point resistances on an infinite-size resistor network can be approximated with those of the 
corresponding finite-size networks. 

In the following section, we give some examples of finite resistor networks such that at the 
limit of the large size of the networks, we obtain some infinite regular networks. Then, we 
approximate the infinite networks with finite ones as illustrated above. 



4 Examples 

In this section we calculate two-point resistances R a p(i), for i — 1,2,3 by using (3-46), (3-55) 
and (3-61), in some important distance-regular networks with diameters d — 1, d — 2 and 
d > 2, respectively. 
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4.1 Complete network K n 

The complete network K n is the simplest example of distance-regular networks. This network 
has n vertices with n(n — l)/2 edges, the degree of each vertex is k = n — 1 also the network 
has diameter d — 1. The intersection array of the network is {c ; 61} = {n — 1; 1}. Clearly, 
this graph has only two strata r (a) = a and i2i(a:) = {/9 : (5 7^ a}. Then, the only two-point 
resistance is R a p(i) which is given by using (3-46) as follows 

2r 1 2r 

R agW = _(l--) = _ for all (3 e iVa). (4-63) 

v — 1 f v 

4.2 Strongly regular networks 

One of the most important distance regular networks are those with diameter d = 2, called 
strongly regular networks. A network with v vertices is strongly regular with parameters 
v, k, A, fj, whenever it is not complete or edgeless and 

(i) each vertex is adjacent to k vertices, 

(ii) for each pair of adjacent vertices there are A vertices adjacent to both, 

(hi) for each pair of non- adjacent vertices there are \i vertices adjacent to both. 
For a strongly regular network, the intersection array is given by 

{c , d; 61, b 2 } = {k, k - A - 1; 1, //}. (4-64) 

One can notice that, if we consider networks with diameter two and maximum degree k 
and a e V, then a has at most k neighbors, and at most k(k — 1) vertices lie at distance two 
from a. Therefore 

v<1 + k + k 2 — k = k 2 + 1, or k> yv — 1, (4-65) 

where, in the following by using the inequality (4-65), we will obtain upper bounds for two- 
point resistances in strongly regular networks. To do so, first we calculate two-point resistances 
for these networks. 



Two-point resistance 19 
By using (3-46), (3-55) and (4-64), we obtain 

Of <?t — i 

= ), and (4-66) 

K V 



2r , . 2k -A 
;(«- A - 1) 

respectively. Then, from (4-65) and (4-66), we obtain the following upper bound for R aj3 (i) 



R *» = K(K -x-i) ( K - > '-— >' (4 - 67) 



< (4-68) 

Now, we consider the following two well-known strongly regular networks. 

A. Petersen network 

Petersen network [13] is a strongly regular network with parameters (v,k, A, rj) = (10,3,0, 1) 
and the intersection array {c , C\, b\, b 2 } = {3, 2; 1,1}. Therefore, by using (4-66) and (4-67), 
we obtain 

3r 4r 
R a/3W = y and R a/3( 2 ) = y- (4-69) 

From (4-69), it is seen that R a po) in Petersen graph saturates the upper bound (4-68). 

B. Normal subgroup scheme 

Definition 2.3 The partition P = {P , Pi, P d } of a finite group G is called a blueprint [13] 
if 

(i) = {e} 

(ii) for i=l,2,...,d if g e Pi then g' 1 e Pi 

(iii) the set of relations Ri = {(a, (3) 6G® G\a~~ l f3 E Pi} on G form an association scheme. 
The set of real conjugacy classes given in Appendix A of Ref. [1] is an example of blueprint on 
G. Also, one can show that in the regular representation, the class sums Pi for i — 0,1, ...,d 
defined as 

Pi= £ 7 eCG, t = 0,l,...,d, (4-70) 
are the adjacency matrices of a blueprint scheme. 
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In [1], it has been shown that, if H be a normal subgroup of G, the following blueprint 
classes 

P = {e}, P 1 = G-{H}, P 2 = H-{e}, (4-71) 

define a strongly regular network with parameters (v, k, A, 77) = (g, g — h, g — 2h,g — h) and 
the following intersection array 

{c , en 61, 62} = - /i, /i - 1; 1, 5 - (4-72) 

where, <? := \G\ and /i := It is interesting to note that in normal subgroup scheme, 

the intersections numbers and other parameters depend only on the cardinalities of the group 
and its normal subgroup. By using (2-13), the QD parameters are given by {a±, a 2 ; loi, ou 2 } = 
{g - 2h,0;g- h, (g - h)(h - 1)}. Also, from (4-71), it is seen that |r 2 (or)| = h - 1. Then, by 
using (4-66) and (4-67), we obtain 

2r(g — 1) 2r 

One should notice that, the maximum degree k for normal subgroup scheme is K max = 9 — 2 
(h = 2), which can be appear in networks with even cardinality such as dihedral group. 
Therefore, for normal subgroup scheme (strongly regular networks with parameters (g, g — 
h,g — 2h, g — h)), we have 

K<g-2, (4-74) 

and therefore, by using (4-65), (4-78) and (4-75) (k = g — h), we obtain upper and lower 
bounds for R al3 w and R a/3 (2) as follows 

Mff - 1) < p < 2ry^T 2r 2r 

^-2) p 5 p 1 

As an example, we consider the dihedral group G = D 2m , where its normal subgroup is 

H = Z m . Therefore, the blueprint classes are given by 

P = { e }, p 1 = {b, ab, a 2 b, a rn ~ l b}, P 2 = {a, a 2 , a (m_1 }, (4-76) 
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which form a strongly regular network with parameters (2m, to, 0, to) and the following inter- 
section array 

{c , ci; 6i, 6 2 } = {to, to - 1; 1, to}. (4-77) 

By using (4-78), we obtain 

tm2?7x — 1 ) 2t 

A*/3(D = 2 ' aIld R <*l#* = ( 4_78 ) 

m p m 

4.3 Cycle network C„ 

A cycle network or cycle is a network that consists of some number of vertices connected in 
a closed chain. The cycle network with v vertices is denoted by C v with k — 2. For odd 
v = 2m + 1, the intersection array is given by 

{c , c m - h, 6 m } = {2, 1, 1; 1, 1, 1}, (4-79) 

where, for even v = 2m, we have 

{c , c m ; h, b m } = {2, 1, 1, 2; 1, 1, 2} (4-80) 

and the network consists of to + 1 strata. We consider v = 2m (the case v = 2m + 1 can be 
considered similarly). Then, by using the recursion relations (2-22) and (2-24) one can obtain 
the following closed form for the Stieltjes function 

G ^ - kWrr ( ^ 81) 

where, T k (x) are Tchebyshev polynomials of the first kind. 
By using (3-46), (3-55) and (3-61) we obtain 

.2m — 1. .to — 1. in ,2m — 3, 

R af3W = r i 2m ), RapW = 2r ^~^~ )' alld Ra(5(V=3r( — ), (4-82) 

respectively. From (4-82), one can easily deduce that 

R af nu)=kr{^^-) k = 1,2,.., m. (4-83) 
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At the limit of the large m, the cycle network tend to the infinite line network and the 
Stieltjes function (4-81) reads as 

G,(x) = -=L=. (4-84) 
\Jx l — 4 

Then, for A we have 

A = lim((x - 2)GJx)) = 0. (4-85) 
Therefore, by using (3-46), (3-55) and (3-61) we obtain 

R a p { i) = r, R af3( 2) = 2r, R a/3W = Sr. (4-86) 

In fact, it can be easily shown that 

R a/3(k) =kr, k = 1,2, ... (4-87) 

where, this result could be obtained from (4-83), for large m. 

As (4-83) indicates, for m larger than ~ 70 the difference \R a pw — kr\ is equal to ~ 0.01, 
where for m larger than ~ 100 we have \R a p(k) —kr\ ~ 0. Therefore, the finite resistor network 
C 2m , with m ~ 100 is a good approximation for the infinite line resistor network. 

4.4 d-cube 

The d-cube, i.e. the hypercube of dimension d, also called Hamming cubes, is a network with 
2 d nodes, each of which can be labeled by an d-bit binary string. Two nodes on the hypercube 
described by bitstrings x and y are are connected by an edge if \x — y\ = 1, where \x\ is the 
Hamming weight of x. In other words, if x and y differ by only a single bit flip, then the two 
corresponding nodes on the graph are connected. Thus, each of the 2 d nodes on the <i-cube 
has degree d. For the ci-cube we have d + 1 strata with 
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The intersection numbers are given by 

bi = d-i, < % < d - 1; a = i, l<i<d. (4-89) 

Then, by using (3-46), (3-55) and (3-61), we obtain 

2 d - 1 2 d - x - 1 

Q/3(1) = ^2^' a/3(2) = (d- l)2 rf -2 r ' 

^ (a) = rf(rf-l)(rf-2) { 2^ } - (4 " 90) 

From (4-90) one can see that, at the limit of the large dimension d, the two-point resistances 
Rapw, i — 1)2,3 tend to zero. Since R a pw tend to zero for d larger than ~ 200. Therefore, 
the finite rf-cube with d larger than ~ 200 is a good approximation for the infinite hypercube 
resistor network. 

4.5 Johnson network 

Let n > 2 and d < n/2. The Johnson network J(n, d) has all d-element subsets of {1, 2, n} 
such that two d-element subsets are adjacent if their intersection has size d—1. Two d-element 
subsets are then at distance % if and only if they have exactly d — i elements in common. The 
Johnson network J(n,d) has v = d ^™-dy. ver ti ces ) diameter d and the valency k = d{n — d). 
Its intersection array is given by 

bi = (d-i)(n-d-i), 0<i<d-l; q = i 2 , l<i<d, (4-91) 

Then, by using (3-46), (3-55) and (3-61), we obtain 

2(n! -d\{n-d)\) 
R ^ W ~ d(n - d)n\ r ' 

2r , , d\(n - d)\(n - 2 - 2d(n - d)) , 

^n/3( 2 ) — ~rr, w ^7 ; r{d(n-d)-(n-2) + -± — ; - -} and 

a/3 d(d- \){n - d)(n - d - l) x v ; v ; n\ 1 

2r 

R ^ (:i> = d(d -l){d-2){n-d){n-d-l){n-d-2) ^ " 2d + 1) + 
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(3w - 2d(n - d) - 10) din - - ^ + [d 2 (n - df - din - d)(3n - 9) - 4(d - l)(w - d - 1) + 

2(n-2)(n-4)](l- rf!(n 7 rf)! )}. (4-92) 

It could be noticed that for a give d, the result (4-92) show that, at the limit of the large 
dimension n, the two-point resistances R a pw, i = 1,2,3 tend to zero. Since R a pm tend to 
zero for n larger than ~ 200. Therefore, the finite Johnson network J(n, d) with n larger than 
~ 200 is a good approximation for the infinite Johnson resistor network. 



5 Two-point resistances in more general networks 

Although, we discussed through the paper only the case of distance-regular networks, but 
also the method can be used for any arbitrary regular network. For calculating two-point 
resistances, we need only to know the Stieltjes function G^(x). For two arbitrary nodes a and 
(3 of the network, we choose one of the nodes, say a, as reference vertex. Then, the Stieltjes 
function G^(x) can be calculated by using the recursion relations (2-22) and (2-24), where, as 
it has be shown in [20], the coefficients and for % — 1, d in the recursion relations are 
obtained by applying the Lanczos algorithm to the adjacency matrix of the network and the 
reference vertex \a). In fact, the adjacency matrix of the network takes a tridiagonal form in 
the orthonormal basis = 0, 1, ...,d} produced by Lanczos algorithm and so we obtain 

again three term recursion relations as (2-22). But, in general, the basis produced by Lanczos 
algorithm do not define a stratification basis, in the sense that, a vertex ket \(3) of the network 
may be appear in more than one of the base vectors \4>i)- In these cases, if d is equal to 
v (the number of vertices of the network), one can write each vertex ket \(5) uniquely as a 
superposition of the base vectors \<f>i) and calculate two-point resistance R a p by calculating 
the entries (0j|L _1 |Q;) for all i — 0,1, d as illustrated through the paper. In the most cases, 
d is less than v. In these cases, we need to obtain some additional orthonormal base vectors 
{\ipi),i = 1, .., v — d — 1} such that the new bases are orthogonal to the subspace spanned 
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by {\<fri),i = 0, l,...,d}. One can obtain some such additional base vectors, by choosing 
a normalized vector orthogonal to the subspace spanned by {\4>i),i = 0, l,...,d} as a new 
reference state and applying the Lanczos algorithm to the adjacency matrix of the network 
and the new reference state. If the number of the new orthonormal base vectors still be less 
than v—d—1, we choose another normalized reference state orthogonal to the subspace spanned 
by all previous orthonormal bases and apply the Lanczos algorithm to the adjacency matrix 
and the new chosen reference state. By repeating this process until to obtain v orthonormal 
basis, one can solve a system of v equations with v unknowns to write each vertex ket (3 as 
superposition of the v orthonormal bases. 



6 Conclusion 

The resistance between two arbitrary nodes in a distance-regular resistor network was obtained 
in terms of the Stieltjes transform of the spectral measure or Stieltjes function associated with 
the network and its derivatives. It was shown that the resistances between a node a and all 
nodes (3 belonging to the same stratum with respect to the a are the same. Also, explicit an- 
alytical formulaes for two-point resistances R a p for (3 belonging to the first, second and third 
stratum with respect to the a were driven in terms of the size of the network and the corre- 
sponding intersection numbers. In particular, the two-point resistances in a strongly regular 
network with parameters (v, k, A, /i) were given in terms of these parameters. Moreover, the 
lower and upper bounds for two-point resistances in strongly regular networks was discussed. It 
was discussed that, the introduced method can be used not only for distance- regular networks, 
but also for any arbitrary regular network by employing the Lanczos algorithm iteratively. 
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Appendix A 

In this appendix, we give the two-point resistances R a p(i), % = 1,2,3 for some important 
distance-regular networks with v < 70. 



The network 


V 


Intersection array 


Ref. 


R a /3W 


R a/ 3(2) 




Icosahedron 


12 


{5, 


2. 


1; 


1. 


2 


5} 


[21] 


llr 
30 


7r 
15 


r 
2 


L (Petersen) 


15 


{4, 


2. 


1: 


1. 


1 


4} 


[21] 


7r 
15 


19r 
30 


2r 
3 


Pappus, 3-cover K 3 3 


18 


{3,2, 


2. 


1; 


1 


1 


2,3} 


[21] 


17r 
27 


8r 
9 


26r 
27 


Desargues 


20 


{3,2,2, 


1. 


1; 


1 


1 


2,2,3} 


[21] 


19r 
30 


8r 
9 


59r 
60 


Dodecahedron 


20 


{3,2,1, 


1. 


1; 


1 


1 


1,2,3} 


[21] 


19r 
30 


9r 
10 


16r 

15 


GH(2, 1) 


21 


{4, 


2. 


2; 


1. 


1 


2} 


[21] 


lOr 
21 


2r 
3 


5r 
7 


Klein 


24 


{7, 


4. 


1: 


1. 


2 


7} 


[21] 


23r 
84 


9r 
28 


r 
3 


GQ(2,4)\ spread 


27 


{8, 


6. 


1; 


1. 


3 


8} 


[21] 


13r 
54 


29r 
108 


5r 
18 


#(3,3) 


27 


{6, 


4. 


2: 


1. 


2 


3} 


[21] 


26r 
81 


31r 
81 


llr 
27 


coxeter 


28 


{3,2, 


2. 


1; 


1 


1 


1,2} 


[21] 


9r 
14 


13r 
14 


73r 
84 


Taylor(P(13)) 


28 


{13, 


6. 


1; 


1 


6 


13} 


[22] 


27r 
182 


44r 
273 


91r 
546 


Tutte's 8-cage 


30 


{3,2, 


2. 


2; 


1. 


1 


1,3} 


[21] 


29r 
45 


14r 
15 


139r 
90 


Taylor (GQ (2, 2)) 


32 


{15, 


8, 


1; 


1 


8 


15} 


[22] 


31r 
240 


llr 
80 


17r 
120 


Taylor(T(6)) 


32 


{15, 


6. 


1; 


1 


6 


15} 


[22] 


31r 
240 


lOlr 
720 


13r 
90 


IG(AG(2A)\pc) 


32 


{4,3, 


3. 


1; 


1 


1 


3,4} 


[23] 


31r 
64 


5r 
8 


125r 
192 


Wells 


32 


{5,4, 


1. 


1: 


1 


1 


4,5} 


[23] 


31r 
80 


15r 
32 


39r 
80 


Hadamard graph 


32 


{8,7, 


4. 


1; 


1 


4 


7,8} 


[24] 


31r 
128 


15r 
56 


249r 
896 


Odd(4) 


35 


{4, 


2. 


1: 


1. 


1 


4} 


[25], [26] 


17r 

35 


22r 
35 


242r 
315 


Sylvester 


36 


{5, 


4. 


4; 


1. 


1 


4} 


[25] 


7r 
18 


17r 
36 


181r 
240 


Taylor(P(17)) 


36 


{17, 


8, 


1; 


1 


8, 


17} 


[22] 


35r 
306 


149r 
1224 


279r 
272 


3-cover Kq^ 


36 


{6,5, 


4. 


1: 


1 


2 


5,6} 


[22] 


35r 
108 


17r 
45 


211r 
540 
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The network 


f 


Intersection array 


Ref. 


R a /3W 


R a /3(2) 




SRG\ spread 


40 


{9, 6,1; 1,2, 9} 


[27] 


13r 
60 


llr 

45 


r 
4 


Ho — Si2(x) 


42 


{16, 5,1; 1,1, 6} 


[25] 


Air 
126 


8r 
21 


34r 
105 


Mathon (Cycl(13,3)) 


42 


{13,8,1; 1,4, 13} 


[28] 


Air 
273 


89r 
546 


91r 
13 


GO{2, 1) 


45 


{4, 2, 2, 2; 1,1, 1,2} 


[23] 


22r 
45 


32r 
45 


Ar 
5 


3-cover GQ(2, 2) 


45 


{6,4,2,1; 1,1,4,6} 


[23] 


44r 
135 


107r 
270 


221r 
540 


Hadamard graph 


48 


{12,11,6,1; 1,6,11, 12} 


[22] 


47r 
288 


23r 
132 


565r 
3168 


IG(AG(2,5)\pc) 


50 


{5, 4, 4,1; 1,1, 4, 5} 


[23] 


49r 
125 


12r 
25 


123r 
250 


Mathon (Cycl(16,3)) 


51 


{16, 10,1; 1,5, 16} 


[28] 


25r 
204 


89r 
680 


68r 
510 


GH(3,1) 


52 


{6, 3, 3; 1,1, 2} 


[29] 


51r 
156 


llr 

26 


23r 
52 


Taylor(SRG(25, 12)) 


52 


{25,12, 1; 1,12, 25} 


[22] 


51r 
650 


319r 
3900 


13r 
156 


3-cover K 9 9 


54 


{9, 8, 6,1; 1,3, 8, 9} 


[22] 


53r 
243 


13r 
54 


239r 
972 


Gosset,Tayl(Schlafli) 


56 


{27,10, 1; 1,10, 27} 


[22] 


55r 
756 


289r 
3780 


49r 
630 


Taylor(Co-Schlafli) 


56 


{27,16, 1; 1,16, 27} 


[22] 


55r 
756 


227r 
3024 


llr 
144 


Perkel 


57 


{6, 5, 2; 1,1, 3} 


[25], [30] 


56r 
171 


22r 
57 


68r 
171 


Mathon(Cycl(ll,5)) 


60 


{11, 8,1; 1,2, 11} 


[28] 


59r 
330 


32r 
165 


r 
5 


Mathon(Cycl(19,3)) 


60 


{19, 12,1; 1,6, 19} 


[28] 


59r 
570 


187r 
1710 


3151r 
570 


Taylor (SRG(29, 14)) 


60 


{29, 14,1; 1,14, 29} 


[22] 


59r 
870 


214r 
3045 


29r 
406 


GH{2,2) 


63 


{6, 4, 4; 1,1, 3} 


[29] 


62r 
189 


76r 
189 


271r 
504 


#(3,4),Doob 


64 


{9, 6, 3; 1,2, 3} 


[22] 


7r 

32 


r 
4 


25r 
96 


Locally Petersen 


65 


{10, 6, 4; 1,2, 5} 




64r 
ozo 


73r 


49r 

lOD 


Doro 


68 


{12, 10, 3; 1,3, 8} 




67r 
408 


145r 
816 


253r 
1020 


Doubled Odd(4) 


70 


{4,3,3,2,2,1,1; 1,1,2,2,3,3,4} 


[22] 


69r 
140 


68r 
105 


869r 
1260 


J(8,4) 


70 


{16,9, 4,1; 1,4,9,16} 


[22] 


69r 
560 


337r 
2520 


691r 
5040 
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